Gravitational wave detection using pulsars: status of the Parkes 
Pulsar Timing Array project 



G. B. Hobbs A , M. Bailes, B , N. D. R. Bhat B , S. Burke- Spolaor*' 3 , D. J. 
Champion A , W. Coles c , A. HotarP , F. Jenet E , L. Kedziora-Chudczer^ , J. 
Khoo A , K. J. Lee E ' F , A. LommerP , R. N. Manchester^ , J. Reynolds^, J. 
Sarkissian A , W. van Straten 3 , S. Tb A,H , J. P. W. Verbiest A,B , D. Yardley^^, 
X. P. You 1 

A Australia Telescope National Facility, CSIRO, P.O. Box 76, Epping, NSW 1710 Australia. 

B Centre for Astrophysics and Supercomputing, Swinburne University of Technology, P.O. Box 218, 

Hawthorn VIC 3122, Australia 

c Electrical and Computer Engineering, University of California at San Diego, La Jolla, California, 
U.SA 

D Department of Imaging and Applied Physics, Curtin University, Bentley, Western Australia, 
Australia 

E Center for Gravitational Wave Astronomy, University of Texas at Brownsville, 80 Fort Brown, 
Brownsville, TX 78520 

F Department of Astronomy, Peking University, 5 Haidian Lu, Beijing, 100871, China 

G Franklin and Marshall College, 415 Harrisburg Pike, Lancaster, PA 17604, U.S.A. 

H Institute of Astronomy, School of Physics A29, The University of Sydney, NSW 2006, Australia 

1 School of Physical Science and Technology, Southwest University, 2 Tiansheng Road, Chongqing 

400715, China 



Abstract: 

The first direct detection of gravitational waves may be made through observations of pulsars. The 
principal aim of pulsar timing array projects being carried out worldwide is to detect ultra-low fre- 
quency gravitational waves (/ ~ 10 -9 -10 -8 Hz). Such waves are expected to be caused by coalescing 
supermassive binary black holes in the cores of merged galaxies. It is also possible that a detectable 
signal could have been produced in the inflationary era or by cosmic strings. In this paper we review 
the current status of the Parkes Pulsar Timing Array project (the only such project in the Southern 
hemisphere) and compare the pulsar timing technique with other forms of gravitational-wave detection 
such as ground- and space-based interferometer systems. 

Keywords: gravitational waves — pulsars: general 



1 Introduction 



Direct detection of gravitational waves (GWs) would 
be of huge importance to the physics and astrophysics 
communities. Detection would simultaneously verify 
our basic understanding of gravitational physics and 
usher in a new era of astronomy. Using models of 
galaxy merg er rates and the properties o f coale scing 
black holes, ISesana. Vecchio fc Colacino I ([2008) and 
references therein predict the existence of an isotropic, 
stochastic, low-frequency GW background. Various 
models of cosmic strings (iDamour fc Vilenkin I l200ll : 
ICaldwell. Battve fc Shellard Ill996l)" and the inflation- 
ary era ijTurner 1119971 : iBovle fc Buonanno 1 12008) also 
predict the existence of a low-frequency GW back- 
ground (|Maggiore I [20001. In the frequency range ac- 
cessible by pulsar timing these backgrounds can be de- 
scribed by their characteristic strain spectrum, h c (f), 



which can be approximated as 

where /i yr = 1/1 yr. The dimensionless amplitude 
of the background, A g , and the spectr al exponent, 
a, d epend on t he typ e of background dJenet et al. I 
120061 ). ISazhin I <| 19781 ) and iDetweilerl <| 1979ft showed 
that these GW signals could be identified through ob- 
servations of pulsars. 

S tandard pulsar timing tech niques (e.g. lEdwards. Hobbs fc Manches 

120061 : fL oruncr fc Kramer II2OO5T ) are used to search for 
the e xistence of GWs. Timing software (Hobbs, Edwards & Manchester 
2006) is used to compare the predictions of a pul- 
sar timing model with measured pulse times of arrival 
(TOAs) . The timing model for each pulsar contains in- 
formation about the pulsar's astrometric, spin and or- 
bital parameters. Any deviations of the actual arrival 



2 



Publications of the Astronomical Society of Australia 



m. 1 

Cosmic 
strings 


1 //$ 


- 






Coalescing binary _ 


Binary 
_ black— holes 
in galaxies 


T Relic \ 
G— waves \ 


/ SN Core 
/ collapse 

1 1 / " 


_ Standard \ 
inflation \ 






Unresolved 
Galactic binaries 


NS-NS ""N/^jj^ 
CoalescenceX^y^Advanced 


PPTA 


LISA LIGO 



-10 



-5 



log 10 (f/Hz) 



Figure 1: Characteristic strain sensitivity for ex- 
isting and proposed GW detectors as a function of 
GW frequency. Predicted signal levels from vari- 
ous astrophysical sources are shown. Figure from 
Hobbs (2008b). 



times from the model predictions - the pulsar timing 
residuals - represent the presence of unmodelled effects 
which will include calibration errors, spin-down irreg- 
ularities and the timing signal caused by GWs. The 
signal from GWs can be disentangled from other un- 
modelled effects by looking for correlations in the tim- 
ing residuals of a large number of pulsars that are dis- 
tributed over the entire sky. Pulsar timing experiments 
are sensitive to GW signals in the ultra-low frequency 
(f ~ 10 -9 -10 -8 Hz) band making them complementary 
to the space- and ground-based interferometers such as 
the Laser Interferometer Space Antenna (LISA) and 
the Laser Interferometer Gravitational Wave Observa- 
tory (LIGO), which are sensitive to higher frequency 
GWs (see Figured] and E}. 

For a given pulsar and GW source, the induced 
timing residuals caused by a GW signal only depend 
on the GW characteristic strain at the pulsar and at 
the Earth. The GW strains at the positions of well- 
separated pulsars will be uncorrelated, whereas the 
component at the Earth will lead to a correlated signal 
in th e timing residuals of all pulsars (Hcllings fc Downs I 
1983). It is this correlated signal that the pulsar tim- 
ing experiments aim to detect. 

The hrst high-precision pulsar timing array exper- 
iments were carried out usin g the Arecibo and Green- 
Bank radio telescopes (e.g. iFoster fc Backer I [l990i . 
Observations of PSRs B1937+21 and B1855+09 were 
subsequently used to o btain a limit on the existence of 
a GW background (e., 



of only 200 ns over a 10 yr period. Timing precisions 
are continuing to impro ve. 

iLommen fc Backer I (|200ll ) described how pulsar 
data sets could be used to constrain the existence of 
binary black holes in the centre of our Galaxy and 
in close-by galaxies. In 2003, a VLBI experiment de- 
tected motions in the radio galaxy 3C66B that were 
thought to b e evidence of a bina r y supermassiy e black 
hole system ijSudou et al. Il2003h . Ijenet et all (|2004h 
determined the expected PSR B1855+09 timing resid- 
uals that would occur due to GW emission from the 
postulated system and compared them with archival 
Arecibo data. The pulsar timing observations allowed 
the existence of the binary system to be ruled out with 
95% confidence. 

The Parkes Pulsar Timing Array (PPTA) project 
began in the year 2003. The three main aims of this 
project are 1) to detect GW signals, 2) to establish a 
pulsar-based time scale for comparison with terrestrial 
time standards and 3) to improve the Solar system 
planetary ephemerides. In this paper, we describe the 
work that has already been carried out towards our 
first aim of GW detection. This includes determin- 
ing the number of pulsars and observing parameters 
required to detect a GW signal, obtaining the most 
stringent constraint on the existence of a GW back- 
ground and producing high quality data sets for use in 
GW astronomy. 



2 Current status of the Parkes 
Pulsar Timing Array project 

The PPTA project has been described in numerous 



papers (for an hi s torical overview see [ Hobbs 



lManchesterll2006l ; iHobbs I l2008bl ; iManchester 



120051 : 



2008). 



In brief, the 20 millisecond pulsars listed in Table [T] 
are observed with an ~ 2-week cadence using the 64- 
m Parkes radio telescope. The pulsars were chosen 
because they 1) only exhibit small-scale timing irreg- 
ularities, 2) have been well studied by earlier projects 
and 3) have a range of separations on the sky. Obser- 
vations are carried out with a 20cm receiver (providing 
256 MHz of bandwidth) and with a 10cm/50cm dual- 
band system (giving 64MHz of bandwidth at 50cm 
and 1GHz of bandwidth at 10cm). Multiple backend 

systems are used including digital filterbanks and the 

Calte ch-Parkes-Swinburne Recorder 2 fCPSR2: lHotan. Bailes fc Ord I 
120061 ) which provides two 64 MHz bands of coherently 
dedispersed data. New backends that are currently be- 
ing commissioned should further improve our timing 
precision. All our o bservations are processed using the 
psrc hive software l|Hotan. van Straten fc Manchester"! 
120041) and pulse TOAs are analysed using TEMPQ2 
Stinebring et al. Il99p| ; iKaspi. Tavl oif HoffibA jpwards fc Manchester 120061 ; lEdwards. Hobbs fc Manchester! 



ll994l ; lMcHugh et al. Hl996l : lLommenll2002h . Unfortu 
nately, the timing precision presented in these early 
papers was too poor to attempt to detect the expected 
GW signals. New pulsar discoveries and better in- 
strumentation has significantl y improved the timing 
precision achievable. Recently IVerbiest et al. I (2008) 
showed that PSR J0437-4715 could be timed using 
the Parkes radio telescope with an rms timing residual 



2006). 



2.1 The initial data sets 

Th e choice of 2 pulsa rs was based on work reported 
bv lJenet et al~l (2005) who determined the minimum 
number of pulsars required to detect the expected GW 
background. This work highlighted the requirement 
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Table 1: Observational parameters for current and future PPTA data sets. N p is the total number of 
observations, T s the data span, a the weighted rms timing residual over the entire data span. For the most 
recent data, ff^cd gi yes the TO A uncertainty, S/N^ d the median S/N and S/N max the maximum S/N. 
S/Nioons is the predicted S/N required to produce a TOA uncertainty of 100 ns and a p is the predicted 
rms timing residual that may be possible with current instrumentation. 
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t The timing residuals for PSR J0437-4715 are currently dominated by systematic effects. Understanding and removing of 
these effects may allow us to reach a TOA timing precision of ~ 20 ns; see text. 



for extremely high precision timing observations - in 
order to detect the GW background within five years 
all 20 puls ars need to have rm s timing residuals of 
~ 100 ns. lYou et al. 1 (|2007aT l showed that the in- 
terstellar medium and Solar wind can produce much 
larger signals in the timing residuals. For instance, 
at an observing frequency close to 1.4 GHz, the So- 
lar wind can produce delays of ~ 100 ns for lines-of- 
sight to the pulsar that pass within 60° of t he Sun and 
much larger delays for closer lines-of-sight (|You et al. I 
2007b]). Fortunately dispersive effects can be removed 
by comparing pulse arrival times at different observ- 
ing frequenci es. An optimal subt raction technique was 
presented bv lYou et al~l (|2007al ). 

Combining the initial PPTA data sets (that s panned 
~2_yr) with archival Arecibo observations allowed Ijenet et 
(2006) to obtain the most stringent limit to date on the 
amplitude of a GW background of A g < 1.1 x 10~ 14 . 
This corresponds to f2 9 [l/(8yr)]/i 2 < 1.9 x 10~ 8 for 
the fractional energy density per unit logarithmic fre- 
quency interval for an astrophysical background (the 
Hubble parameter H = lOO/i kms' 1 Mpc" 1 ). This 
result constrained the merger rate of supermassive bi- 
nary black-hole systems, ruled out some relationships 



between the black-hole mass and the galactic halo mass, 
constrained the rate of expansion in the inflationary 
era and provided an upper bound on the dimensionless 
tension of a cos mic-string backgroun d. Unfortunately, 
the technique of I Jenet et al. I j|2006T ) is only valid for 
data sets that have a "white spectrum" (defined as 
being a spectrum that is independent of frequency). 
The initial data sets, and previously published work 
fe.g.lKa spi. Tavlor fc R vballl 994: Cog nard fc Backer I 
l2004l ; ISplaver et al. l2005l ; IVerbiest et al. 120081 ') . clearly 
show that many millisecond pulsar timing residuals 
are affected by unexplained timing irregularities which 
may be related to the timing noise phenomena s een in 
younger pulsars (jHobbs. Lvne fc Kramer I [2006) . un- 
modelled instrumental eff ects, poor remova l of inter- 
alslj ellar medium effects fe.g lStinebring |[2006l l or intrin- 
sic pulse shape variability. This limitation restricted 
our initial analysis to only seven of the 20 PPTA pul- 
sars. We have recently developed a new technique that 
makes no assumption on the spectrum of the timing 
residuals. This work will be presented in a forthcom- 
ing paper. 
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3 Towards detecting GW sig- 
nals 

Timing residuals for our current PPTA data sets are 
summarised in Table [1] To form these data sets we 
have selected the backend systems providing most pre- 
cise TOAs. The earliest data are usually from the 
CPSR2 coherent dedispersion system and the most re- 
cent data from the digital filterbank systems. Columns 
2, 3 and 4 provide the total number of observations, 
the time span and the weighted rms timing residual 
over our entire data span respectively (obtained after 
fitting only for the pulsar's astrometric and orbital pa- 
rameters along with the pulse frequency and its first 
derivative). These timing residuals have not been cor- 
rected for the effects of dispersion measure variations 
caused by the interstellar medium. Note also that mul- 
tiple, independent observations at different frequencies 
have been made during the same observing session. A 
weighted average (in 2-weekly intervals) of the resid- 
uals for PSR J0613-0200 gives an rms of 842 ns with 
102 data points (cf. 1.1 (JB listed in the table). As 
our observing systems have been improved over this 
period, the median arrival-time uncertainties for our 
most recent observations are listed in Column 5 of Ta- 
ble ED 

Even though many of these data sets have lower 
rms values than have been achieved in the past, we 
are still far from the goal of ~ 100 ns for the majority 
of our sample. For a given backend system it is possible 
to use a noise-free approximation of the pulse profile 
to determine the best timing precision achievable for a 
given pulsar. This can be carried out thr ough simula- 
tions or by using analytical ex pressions l|van Straten I 
l2ood : Id owns fc Reichlev II1983T ). A typical example is 
shown in Figure [2] where the actual arrival time uncer- 
tainties for PSR J1713+0747 are compared with the 
theoretical prediction. For this pulsar it is clear that 
we should be able to achieve typical arrival time uncer- 
tainties ~ 100 ns. PSR J0437-4715 provides our best 
quality data sets. However, as shown in Figure [2] we 
are far from the theoretical timing precision achievable 
of ~ 20 ns for this pulsar because of un-explained sys- 
tematic effects. In Table [1] we list the signal-to-noise 
(S/N) ratio of an individual observation that would be 
required in order to measure an arrival time with a 
precision of 100 ns (S/Nioo ns ). The table also contains 
the median (S/N^d) an( i maximum (S/N™^) S/N ra- 
tios achieved to date using our most recent observ- 
ing system. Further improvements should be possible 
by improved calibration, using a new coherent dedis- 
persion system currently being commissioned, longer 
observation lengths and by obtaining pulse times-of- 
arrival using the information av ailable from the Stokes 
parameters (jvan Straten 1120061 ) instead of considering 
only the total intensity profile. 

In order to achieve rms timing residuals of ~100 ns 
over five years of observations it is necessary that the 
pulsars are intrinsically stable over this time scale. A 
generalisation of the Allan variance h as been developed 
for an alysing pulsar timing residuals (iMatsakis. Taylor fe 
1997). This measure of stability as a function of time- 
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Signal/ noise 

Figure 2: Arrival time precision versus signal to 
noise ratio for PSRs J1713+0747 and J0437-4715. 
For PSR J1713+0747 the measured values are 
plotted as circles and a theoretical estimate is 
drawn as the solid line. For PSR J0437-4715 we 
use 'cross' symbols and a dashed line respectively. 
The observed weighted rms timing residuals are 
indicated by dashed lines. The dotted lines repre- 
sent timing precisions of 100 ns and 500 ns. 

span, u z (t), has been plotted in Figure [3] for our best- 
timed pulsar, PSR J0437— 4715, and for one of our 
poorer pulsars, PSR J1045— 4509. The figure also in- 
cludes measurements of a z for the difference between 
the realisations of terrestrial time published by the Na- 
tional Institute of Standards and Technology (NIST) 
and the Physikalisch-Technische Bundesanstalt (PTB). 
These initial results demonstrate that our best pulsars 
are stable enough to search for irregularities in the ter- 
restrial time scales and could be used to form a pulsar- 
based time scale. 

3.1 The GW background 

A GW background will induce timing residuals that 
are correlated between different pulsars. The general 
relativistic prediction for t his correlation is shown a s 
the solid line in Figure H ([Heflings fc Downs I I"l983l ). 
Determining the correlation coefficients for our data 
sets is not trivial due to the uneven sampling, differing 
data lengths and the presence of low-frequency noise. 
However, an initial analysis in which we averaged the 
measured timing residuals into 2-weekly sessions be- 
fore using standard correlation techniques, produced 
the correlation coefficients shown in Figure Clearly 
we have not yet made a significant detection. 

A reference tim ing array project was defined by 
IJenet et al. I (|2005h containing 20 pulsars, timed with 
an rms timing residual of 100 ns with weekly sam- 
pling over 5 years. Such observations would give a 
detection significance of ~ 3 for a background with 
A g = 2 x 10 -15 . With our current rms timing resid- 
uals an d fortnightly sampling we would require data 
EuijsHiksi ig more than 16 yr in order to achieve this sen- 
sitivity. For rough estimates of the timing precision 
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Figure 3: The <j z (t) stability parameter for PSRs 
J0437-4715 (solid circle) and J1045-4509 (star 
symbols). For comparison the difference between 
two realisations of the terrestrial time scale is also 
shown (open circles). Error bars are only plotted 
if they are larger than the size of the data point. 



Figure 4: Correlation between all pairs of pulsars 
plotted versus their angular separation. The solid 
line is the predicted correlation curve due to a GW 
background. 



that we are likely to achieve with improved instrumen- 
tation (listed in column 9 in Table [T] and obtained by 
postulating possible future S/N ratios and their corre- 
sponding TO A errors), we require data spans of ~10 yr 
to reach this sensitivity. 

In order to reduce the time required to create such 
data sets, we are developin g collaborations with the 
Nort h AmericarQ E uropean ( Janssen et al. 2008; Sta ppers 
l2006ft and Chinese ijWang et al. 1 120021 ) timing array 
projects. It is also expected that on-going pulsar sur- 
veys will discover new, stable pulsars that can be used 
in future timing array projects. In the longer term we 
expect that the Square Kilomet r e Array (SKA) tele- 
scope (ICarilli fc Rawlingsl |2004| ; ICordes et aj~l |2004 : 
iKramer 112001 7 will provide high quality data sets for 
many hundreds of pulsars. Such observations will al- 
low the properties of the GW background to be studied 
in detail. For instance, predictions of various theories 
of gravity for the polarisation properties of GWs could 
be tested (Lee et al., in press). 



binary black-hole systems as a function of black-hole 
mass and redshift. 

Sources of burst GW emission that may be de- 
tectable with current and proposed pulsar timing ex- 
periments i nclude 1) the formation of supermassive 
black holes ijThorne fe Braginskii I I 1976T ) , 2) highly ec- 
centr ic supermassive black hole binaries llEnoki fc N agashi ma | 
20071 ') , 3) close encounters of mass ive objects jKoc^isl^aspar" fc Markal 
2006) and 4) cosmic string cusps ijDamour fc Vilenkin I 
2DuLi)-"j Our ability to detect and localise the position of 



3.2 Single GW sources 

Single GW sources can be divided into individual su- 
permassive black-hole binary systems and burst GW 
sources. Order-of-magnitude calculations demonstrate 
that a system with a chirp mass of 10 9 Mq and an or- 
bital period of 10 yr at a distance of 20 Mpc will in- 
duce a sinusoidal signal in the timing residuals with 
amplitude ~ 100 ns at a frequency of twice the orbital 
frequency of the binary black holes. More detailed 
calculations, which take various cosmological effects 
into account when determining our sensitivity to sin- 
gle GW sources, will be published shortly. This new 
work demonstrates that our data sets can be used to 
place constraints on the merger rate of supermassive 



burst sources depends on our having accurate data sets 
for a large number of pulsars. By sharing data with 
timing array projects in the Northern hemisphere we 
expect to improve both the angular resolution of our 
detector and our sensitivity to burst sources. 

4 The pulsar timing GW band 

Any GW signal with a period longer than the time- 
span of the data is largely absorbed by the fitting pro- 
cess. Pulsar timing is therefore only sensitive to GW 
signals with periods less than the time-span of the 
data (/ < 10~ 8 Hz). The planned space-based GW 
detector, LISA, is sensitive to higher frequency GWs 
(/ ~ 10 -3 Hz) and the ground based interferometers 
to even higher frequencies (/ ~ 10 2 Hz). In general, 
the expected GW sources are different for each type of 
detector. However, some models of cosmic strings and 
the inflationary era produce GW backgrounds which 
could, pote ntially, be detect able at all GW observing 
frequencies (iMaggior e l|2000|) 



ihttp:/ /nanograv.org 



As noted by Sesana. Vecchio fc Colacino I (2008), 
detection of a background from coalescing binary black 
hole systems by both the pulsar timing array experi- 
ments and LISA would provide strong constraints on 
models of the ass embly of mass i ve bi nary black hole 
systems. Recentlv. lPitkin et al. I l|2008f ) discussed other 
areas where pulsar timing and LISA may provide com- 
plementary observations. For instance, they consider 
the possibilities of studying the inspiral stage of a merger 
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using pulsar observations and the ring-down stage us- 
ing LISA. 



5 Conclusions 

Current pulsar timing array experiments have the po- 
tential to detect low-frequency GW signals within a 
decade. Future experiments with the Square Kilome- 
tre Array telescope should allow detailed observations 
of both a stochastic GW background and individual 
GW sources. The pulsar timing technique is comple- 
mentary to other projects that are attempting to de- 
tect much higher frequency GWs. 
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